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Abstract 

We present results of weak leasing cluster counts obtained from 11 degree^ Subaru/SuprimeCam data. 
Although the area is much smaller than previous work dealing with weak leasing peak statistics, the number 
density of galaxies usable for weak leasing analysis is about twice as large as those. The higher galaxy 
number density reduces the noise in the weak leasing mass maps, and thus increases the signal-to-noise 
ratio {S/N) of peaks of the leasing signal due to massive clusters. This enables us to construct a weak 
leasing selected cluster sample by adopting a high threshold S/N, such that the contamination rate due 
to false signals is small. We find 6 peaks with S/N > 5. For all the peaks, previously identified clusters 
of galaxies are matched within a separation of 1 arcmin, demonstrating good correspondence between the 
peaks and clusters of galaxies. We evaluate the statistical error in the weak leasing cluster counts using 
mock weak leasing data generated from full-sky ray-tracing simulations, and find Apeak = 6±3.1 in an 
effective area of 9.0 degree^. We compare the measured weak leasing cluster counts with the theoretical 
model prediction based on halo models and place the constraint on Ulm — erg plane which is found to be 
consistent with currently standard ACDM models. It is demonstrated that the weak leasing cluster counts 
can place a unique constraint on erg — cq plane, where cq is the normalization of the dark matter halo 
mass-concentration relationship. Finally we discuss prospects for ongoing/future wide field optical galaxy 
surveys. 

Key words: cosmology: observations — dark matter — cosmological parameters — large-scale struc¬ 
ture of universe 


1. Introduction 

The A-flat universe with cold dark matter (ACDM) is 
now considered as the standard theoretical framework for 
cosmic structure formation. In the ACDM paradigm, dark 
matter haloes of a galaxy and a cluster of galaxies form 
though the hierarchical bottom-up assembly of smaller 
structures. The averaged density profile of the dark mat¬ 
ter haloes found in A-body simulations has a universal 
broken power-law form (Navarro et al. 1997). The univer¬ 
sal profile is characterized by a single parameter (called 
the concentration parameter) describing the ratio of the 
scale radius (where the power-law slope is approximately 
—2) to the virial radius which defines a halo’s mass. The 
concentration parameter is known to be related to the 
halo mass, a correlation which most likely arises from the 
mass assembly history of haloes (Bullock et al. 2001; Zhao 
et al. 2009; Ludlow et al. 2014; Diemer & Kravtsov 2015). 
Since the mass assembly history depends on the cosmolog¬ 
ical model, the mass-concentration relation does as well 
(Maccio et al. 2008; Kwan et al. 2013). Therefore plac¬ 
ing constraints on both the cosmological model and the 
mass-concentration relation simultaneously allows us to 
test the ACDM structure formation scenario. 

Clusters of galaxies serve as one of the most important 
tests of the cosmic structure formation model. Due to 
their high virial temperatures, the dissipative cooling of 
baryons is less effective, and thus the density profile of 


clusters may be well approximated by the dark matter 
distribution, allowing us to test the mass-concentration 
relation (see for a recent review Ettori et al. 2013, and 
references therein). Also, cluster number counts have 
placed useful constraint on the cosmological parameters 
(Vikhlinin et al. 2009; Mantz et al. 2010; Rozo et al. 2010). 

It is argued in recent studies that number counts of 
clusters of galaxies identified in weak leasing mass maps 
(which we call weak leasing cluster counts) can be a pow¬ 
erful probe of both the cosmological model and the mass- 
concentration relation (Mainini & Romano 2014; Cardone 
et al. 2014). This can be simply understood as follows: 
First, the number density of massive clusters of galax¬ 
ies, which appear as high peaks on the weak leasing 
mass maps, depends of the cosmological model (Borgani 
et al. 2001; Henry 2004). Second, peak heights depend 
on the density profile of clusters and thus on the mass- 
concentration relation (King & Mead 2011; Hamana et al. 
2012 ). 

Both observational techniques and theoretical models of 
weak leasing cluster counts have been greatly developed 
in the last decade. On the observational side, system¬ 
atic weak leasing cluster searches have become practicable 
(Miyazaki et al. 2002a, 2007; Wittman et al. 2006; Gavazzi 
& Soucail 2007; Schirmer et al. 2007; Kubo et al. 2009; 
Bellagamba et al. 2011; Shan et al. 2012). Miyazaki et al. 
(2007) presented the first large sample of weak leasing 
peaks (which are plausible candidates of clusters/groups 
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of galaxies) in 16.7 degree^ Subaru weak lensing survey 
data, with 100 peaks with signal-to-noise ratio exceeding 
S/N > 3.7. Shan et al. (2012) reported 301 weak lens¬ 
ing high peaks {S/N > 3.5) located from 64 degree^ of 
CFHTLS-wide data, among those peaks, they confirmed 
85 groups/clusters. On the theoretical side, models of the 
weak lensing cluster counts have been developed: Based 
on the pioneering work by Bartelmann et al. (2001), some 
improvements have been made, including the effect of 
the galaxy intrinsic ellipticities (Hamana et al. 2004; Fan 
et al. 2010), the effect of the large-scale structures (Marian 
et al. 2010; Maturi et al. 2010), the effect of the diversity 
of the dark matter distribution within clusters (Hamana 
et al. 2012), and an optimal choice of the filter function 
(Hennawi & Spergel 2005; Maturi et al. 2005). Now, the 
theoretical models were tested (or calibrated) against N- 
body simulations and good agreements were found (see 
the above references). 

In fact, weak lensing peak counts have now become 
a practical tool for investigating cosmological models. 
Liu et al. (2014a) used weak lensing peak counts com¬ 
bined with weak lensing power spectrum measured from 
the CFHTLenS data to constrain cosmological parame¬ 
ters. They used not only high peaks but also low peaks 
(—0.04 < K < 0.12) in order to efficiently utilize cosmo¬ 
logical information contained in the weak lensing mass 
map. They found that the constraints from peak counts 
are comparable to those from the cosmic shear power 
spectrum. Liu et al. (2014b) used weak lensing peak 
counts within 3 < S/N < 6 from ~ 130 degree^ of the 
Canada-France-Hawaii Telescope Stripe82 Survey data. 
They demonstrated the capability to put a constraint on 
the mass-concentration relation together with the cosmo¬ 
logical parameters from weak lensing peak counts. 

In this paper, we examine weak lensing cluster counts 
obtained from 11 degree^ of weak lensing mass maps gen¬ 
erated from deep Subaru/SuprimeCam data. In contrast 
to previous work (Liu et al. 2014a,b), we deal only with 
high peaks, with S/N > 5. Such high peaks mostly asso¬ 
ciate with a real massive cluster of galaxies, and thus we 
shall use the term “weak lensing cluster counts” instead 
of “weak lensing peak counts”. Our higher threshold en¬ 
ables us to make a reliable theoretical prediction of the 
weak lensing cluster counts utilizing models of dark mat¬ 
ter haloes which are well calibrated with iV-body simula¬ 
tions. For the first time, we use the weak lensing cluster 
counts to place cosmological parameters and the mass- 
concentration relationship. In doing this, we explore, in 
an experimental manner, the capability and current issues 
of weak lensing cluster counts for cosmological studies. 

The structure of this paper is as follows. In Section 2, 
we summarize the Subaru/SuprimeCam data, data pro¬ 
cessing, and weak lensing shape measurement. We also 
examine the additive and multiplicative bias in the mea¬ 
sured galaxy ellipticities. In section 3, we first describe 
construction of weak lensing mass maps. Then we exam¬ 
ine statistical properties of peaks paying special attention 
to the residual systematics associated with data reduction 
and/or weak lensing shear measurement. In that section, 
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Table 1. Summary of data sets. See Fig. 4 for the coordinate 
and shape of each field 


Field name 

Area“ 

rig'’ 

Seeing'’ 

XMM-LSS 

3.6 (2.8) 

23.7 (20.9) 

0.51-0.71 

COSMOS 

2.1 (1.6) 

28.8 (26.3) 

0.50-0.53 

Lockman-hole 

2.1 (1.6) 

25.5 (23.8) 

0.42-0.55 

ELAIS N1 

3.6 (2.9) 

24.6 (22.0) 

0.45-0.69 


“ Area after masking regions affected by bright stars in 
unit of degree^. The numbers in parentheses are the ef¬ 
fective area (A®-^-^) after cutting the edge regions within 
1.5 arcmin from the boundary. 

^ The mean number density of galaxies used for the weak 
lensing analysis in the unit of arcmin“^. The numbers in 
parentheses are the effective number density defined by 

(Heymans et al. 2012). 

Seeing full width at half-maximum in the unit of arcsec. 
The smallest and largest values among the pointings are 
given. 

we also summarize the high S/N peaks and their pos¬ 
sible counterpart clusters of galaxies searched in known 
cluster databases. In Section 4, the measured weak lens¬ 
ing cluster counts are compared with the theoretical pre¬ 
diction to place constraints on the cosmological param¬ 
eters and the dark matter halo mass-concentration rela¬ 
tionship. Finally, summary and discussion are given in 
Section 5. In Appendix 1, the method for conversion of 
WCSs from the astrometric SIP to TPV convention, which 
is used in the present work, is described. 

2. Data analysis 

2.1. Data set and data processing 

We collected F-band data taken with the 
Subaru/SuprimeCam (Miyazaki et al. 2002b) from 
the data archive SMOKA^ (Baba et al. 2002), under the 
following three conditions: data, which are composed of a 
number of pointings, are contiguous at least 2 degree^, the 
exposure time for each pointing is longer than 1800 sec, 
and the seeing full width at half-maximum (FWHM) of 
each exposure is better than 0.7 arcsec. Four data sets 
meet these requirements and are summarized in Table 1. 
Data of each pointing are composed of several dithered 
exposures. Dithering scales are a few arcmin except for 
the COSMOS field for which a specific dithering pattern 
was adopted (a combination of dithered pointings of the 
normal “north is up” exposures and 90 degree rotated 
exposures, see Taniguchi et al. 2007). 

Each CCD dataset was reduced using an image anal¬ 
ysis pipeline, hscPipe, developed for Hyper SuprimeCam 
(Miyazaki et al. 2012) data (Furusawa et al. in prep.), 
which is based on LSST-stack (Ivezic et al. 2008; Axelrod 
et al. 2010), with being tuned for SuprimeCam data. The 
mosaic stacking was done with hscPipe using the standard 
correlation technique between star positions on CCD co¬ 
ordinates and an external star catalog (we used the SDSS 

^ http://snioka.nao.ac.jp/ 
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Fig. 1. The number densities of galaxies used for weak lens¬ 
ing analysis are plotted against the seeing FWHM for each 
pointing. The different symbols represent the different expo¬ 
sure time. 


DR8), and also between the star positions among differ¬ 
ent exposures. We made stacked data on a pointing-to- 
pointing basis except for COSMOS data for which a single 
stacked data was made by using SCAMP (Bertin 2006) 
and SWarp (Bertin et al. 2002) from all the CCD data. 
Object detections were performed with SExtractor (Bertin 
& Arnouts 1996). Stars are selected in a standard way 
by identifying the appropriate branch in the magnitude- 
object size (FWHM) plane, along with the detection sig¬ 
nificance cut S/N > 10. The number density of stars is 
found to be ~ 1 arcmin”^ for the four helds. We only 
use galaxies met the following three conditions, (i) the 
detection significance of S/N > 8, (ii) FWHM is larger 
than the stellar branch, and (iii) the AB magnitude (for 
which we adopted SExtractor’s MAG_AUT0) is in the range 
of 22 < i' < 25. The number density of resulting galaxy 
catalogue varies among pointings due mainly to the vari¬ 
ation in the seeing and weakly to the exposure time (see 
Fig. 1). 

2.2. Weak lensing shear estimate 

For the weak lensing shear estimates, we employ lensfit 
(Miller et al. 2007; Hitching et al. 2008; Miller et al. 2013), 
which is a model-fitting method. We refer the reader to 
the above references for details of its method, implemen¬ 
tation, and results from tests on image simulations. One 
important function of lensfit which should be noticed is 
that galaxies on each individual exposure image are fit¬ 
ted by models that have been convolved with the PSF for 
the image, and the resulting likelihoods multiply to obtain 
the final combined likelihood. This can avoid problematic 
issues in measuring weak lensing shears from a stacked 
image; including deformation of object shapes by interpo¬ 
lation of pixel data needed for resampling onto a co-added 


pixel grid (Rhodes et al. 2007; Hamana & Miyazaki 2008; 
Miller et al. 2013), and a complex PSF variation over a 
stacked image caused by co-addition of different individual 
PSFs (Hamana et al. 2013). 

Application of lensfit to SuprimeCam data was first 
done by Dietrich et al. (2012). We used the software 

suite which is designed for dealing with SuprimeCam data. 
For the actual implementation of lensfit, we basically fol¬ 
low one developed for CFHTLenS (Miller et al. 2013) and 
thus the method is likelihood-based, and we adopt the 
same priors for the galaxy model. In order to adjust 
lensfit to SuprimeCam data, we have made following two 
modifications. (1) The size of postage-stamp image for 
galaxies was set to be 40 x 40 (48 x 48 for CFHTLenS) 
because of the following two reasons: The first is the 
difference in the CCD pixel sizes between SuprimeCam 
(0.202 arcsec) and MegaCam on CFHT (0.185 arcsec), 
and second is that more smaller galaxies are detected in 
SuprimeCam data owing to their deeper image depth. 
(2) The definition of a boundary of objects was modi¬ 
fied in the following way. In the analysis of CFHTLenS 
(Miller et al. 2013), 2cr isophote was adopted for a bound¬ 
ary of objects in the following two processes; (I) masking 
neighbour objects in the same postage-stamp image, and 
(H) clipping of possible blending objects identified by the 
overlapping of isophotes of close objects and by a mis¬ 
match between the peak and centroid positions. Since the 
noise level of an image depends on the depth of data and 
our data are much deeper than CFHTLenS data (thus 
smaller a), adopting the same 2a isophote extends the 
object boundary to outer skirt and significantly increases 
a chance a object to be classified as a blending object 
and to be rejected. We thus decided to set the threshold 
isophote (/thO") being the CFHTLenS-2(7 corresponding 
value; /th = 2^Texp/^QQsec, where Texp is the total expo¬ 
sure time for each pointing and 900sec is the SuprimeCam 
corresponding exposure time of the CFHTLenS (3600sec). 
Galaxies that are excluded in the process of lensfit shape 
measurement is mostly due to the above criteria, plus a 
small fraction of failures in the model-fitting step. The 
fractions of excluded galaxies vary among pointings, and 
the average value is 0.22. This is comparable to that for 
the CFHTLenS (^ 0.2), supporting that our exposure¬ 
time dependent threshold results in a similar level of 
screening of large, blended or complex galaxies. 

In lensfit, the point spread functions (PSFs) are rep¬ 
resented as postage-stamp data of pixel values measured 
from images of stars. The PSFs are measured for each 
exposure. The spatial variation of PSFs (pixel values of 
PSF postage-stamps) are modeled by a two-dimensional 
polynomial function of position in the SuprimeCam’s field- 
of-view (FoV) that consists of a mosaic of 2 x 5 CCD chips 
with low order coefficients being allowed to vary between 
CCD chips. We employ a third-order polynomial for the 
FoV-based model and a second-order for the chip-based 
model. We did not find meaningful improvement using 
fourth-order for the FoV model. 

We adopted the same galaxy weighting scheme as one 
defined in Miller et al. (2013), which is according to the 
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Fig. 2. The distribution of lensing weight for 8549 measured 
galaxies in one typical pointing, ELAIS_02 (seeing FWHM 
~ 0.57 arcsec and Texp = 3600 sec). The weights are shown as 
a function of (left) i'-band magnitude, (center) fitted semi¬ 
major axis disk scalelength, and (right) signal-to-noise ratio. 


width of the likelihood surface in model-fitting parameter 
space. Fig. 2 shows an example of distribution of lens¬ 
ing weights in one typical pointing, as a function of i'- 
band magnitude, fitted semi-major axis disk scalelength, 
and signal-to-noise ratio. Comparing with the correspond¬ 
ing plots for CFHTLenS (Fig. 2 of Miller et al. (2013)), 
one may find very similar trends between two analyses, 
confirming that the weighting scheme properly works on 
SuprimeCam data. 

It is known that galaxy shears estimated by lensfit 
are affected by multiplicative and additive bias, which 
are usually represented by (Heymans et al. 2012), e°'^® = 
(1-I-m)[ 7 -h e™*]-I-c. Miller et al. (2013) used realistic 
image simulations to estimate the multiplicative bias in 
CFHTLenS data, and found a size- and 5'/A^-dependent 
bias with an average value of 6 percent. They derived the 
fitting function of the multiplicative bias [equation (14) of 
Miller et al. (2013)] that we employ for the current study. 
Using the fitting function, we evaluated an weighted aver¬ 
age bias over source galaxies and found (m) = —0.059. It is 
not clear whether the fitting function that was derived us¬ 
ing mock CFHTLenS simulation data is applicable to the 
SuprimeCam data. We may, however, say that this does 
not seriously affect our results for the following two rea¬ 
sons: First, we don’t use a lensing induced quantity (shear 
or convergence) directly, but use the signal-to-noise ratio 
map of the lensing convergence (the convergence map nor¬ 
malized by its noise map). Therefore, our primary obser¬ 
vational quantity, the peak counts, is not affected by the 
mutiplicative bias (see the next section). We, however, 
need to take the multiplicative bias into account in mak¬ 
ing the theoretical prediction of the peak counts. In order 
to do this, we will make a crude bias correction utilizing 
the average value of (m) = —0.059 (see §4.1 for details). 
Since the statistical error in the peak counts is estimated 
as about 50 percent (see §4.2), the uncertainty associated 
with the multiplicative bias should be minor, that is the 
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second reason. 

Concerning the additive bias, we evaluated it in an 
empirical manner using weighted mean value of ellip- 
ticities over galaxies and found non-zero values with; 
(ci) = (tcef") = -0.00046±0.00030, and (ca) = (tcef") = 
0.0012 ± 0.00030. In order to look into details of this, we 
computed the additive bias as a function of i'-band mag¬ 
nitude and scalelength and found that the non-zero bi¬ 
ases are mainly caused by faint-small galaxies as is shown 
in Fig. 3. This is in marked contrast to the result of 
CFHTLenS that the small-bright galaxies were found to 
be a cause of the significant additive bias in ( 02 ) (Heymans 
et al. 2012 , note that their (ci) = 0 . 0001 ± 0.0001 was found 
to be consistent with zero). The cause of these non-zero 
bias is unknown, but those different trends for different 
data sets may suggest that it may be caused by any instru¬ 
mental or data-dependent effect. Heymans et al. (2012) 
calibrated the bias in an empirical manner, but we don’t 
do it because of the following two reasons: The first is that 
we search for high peaks in weak lensing mass maps which 
are obtained from azimuthally averaged tangential galaxy 
shears with an axisymmetric kernel. In that operation, the 
additive bias is partly canceled out. The second is that we 
deal with very high peaks corresponding to 7 ^ 0.08, thus 
the additive bias with the above small amplitude does not 
have a significant impact on our results. 


3. Weak lensing mass maps and high S/N peaks 

3.1. Weak lensing mass maps 

The weak lensing mass map which is the smoothed lens¬ 
ing convergence field (k) is evaluated from the tangential 
shear data by (Schneider 1996) 

ic{e) = (i) 


where 7 t( 0 : 6) is the tangential component of the shear at 
position (p relative to the point 6 and Q is the filter func¬ 
tion for which we adopt the truncated Gaussian function 
(for K field) (Hamana et al. 2012), 


Q{e) 



( 2 ) 


for 9 < 9o and Q = 0 elsewhere. A nice property of this 
hlter that should be noticed is that it has less power on the 
inner region {9 < 9c, see Fig. 1 of Hamana et al. (2012)) 
thus a Af-peak resulting from a cluster of galaxies is less 
affected by contamination from cluster member galaxies. 

In an actual computation, /C is evaluated on regular grid 
points with a grid spacing of 0.15 arcmin using equation 
( 1 ), but the integral in that equation is replaced with the 
summation over galaxies taking into account the lensing 
weight. The root-mean-square (RMS) noise coming from 
intrinsic ellipticity of galaxies (which we call the galaxy 
shape noise and denote by Ushape) is evaluated on each 
grid point in the following way: A “noise K, map” is eval¬ 
uated in the same way as the mass map but the orienta¬ 
tion of galaxies is randomized. 1000 realizations of “noise 
/C maps” are produced, and the RMS of 1000 /C values 
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Fig. 3. Top-panels: Mean ellipticities of galaxies as a function of semi-major axis desk scalelengtli in three f'-band magnitude 
ranges. Open circles and triangles are for the ei and 62 component of the galaxy ellipticities, respectively, obtained from lensfit 
model-fitting. Bottom-panels: The weighted galaxy fraction for each magnitude-scalelength grid. 


is evaluated on each grid, which we take as the local es¬ 
timate of CTshape- The signal-to-noise ratio of weak leas¬ 
ing mass map is defined by S/N = /C/o-ghape- The filter 
parameters should be chosen so that the S/N is maxi¬ 
mized for expected target clusters (i.e. M > 
at 0.1 < z < 0.5, see Hamana et al. (2004)). We take 
Og = 1-5 arcmin and 9o = 15 arcmin. We show in Fig. 4 
the IC-S/N maps for four fields. Note that a K.-S/N map is 
unaffected by the shear multiplicative bias (m) because we 
used the same galaxy ellipticity data for evaluating both 
the signal and noise, and thus the effect of the multiplica¬ 
tive bias is canceled out. Since in the following analysis 
we use only S/N values, we did not correct for the shear 
multiplicative bias. 

On and around regions where no source galaxy is avail¬ 
able due to imaging data being affected by bright stars 
or large nearby galaxies, K, may not be accurately eval¬ 
uated. We define “data-region”, “masked-region” and 
“edge-region” by using the distribution of source galax¬ 
ies as follows: First, for each grid, we check if there is a 
galaxy within 0.6 arcmin from the grid center. If there is 
no galaxy, then the grid is flagged as “no-galaxy”. After 
performing the procedure for all the grids, all the “no¬ 
galaxy” grids plus all the grids within 0.6 arcmin from all 
the “no-galaxy” grids are defined as the “masked-region”. 
All the grids located within 1.5 arcmin (we take this value 
by setting it equal to 9g) from any of “masked-region” 
grids are defined as the “edge-region”. All the rest of 
grids are defined as the “data-region”. Areas of four fields 
are summarized in Table 1. The total survey area (“data- 
region”-!- “edge-region”) is 11.4 degree^ with the “edge- 
region” accounting for 21 percent. 

Before looking into very high peaks, we examine statis¬ 
tical properties of peaks paying a special attention to an 
impact of residual systematics associated with data reduc- 
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Fig. 5. The probability distribution function of positive- 
and negative-curvature peaks (peaks and troughs) is shown. 
Top-panel is for K map, whereas bottom-panel is for the 
B-mode K. map. The histogram shows the results from ob¬ 
served K and ICb maps. The points with error bars show the 
mean and RMS of PDF measured from 1000 random realiza¬ 
tions of “noise AC maps”, which represents the PDF expected 
for pure shape noise data. 
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Fig. 4. K. maps for the XMM-LSS (top), COSMOS {middle-left), Lockman-hole {middle-right), and ELAIS-Nl {bottom) are shown 
by blue scale. Gray regions are either outside of data regions or masked regions where the data are affected by bright stars. High 
peaks with S/N > 4.5 are marked with red circle along with ID given in Table 2. Black circles are for high peaks within the edge 
region. 
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No. ] 

tion and/or lensing shear measurement. We computed the 
probability distribution function of positive- and negative- 
curvature peaks (peaks and troughs) from the 1C maps, 
and plotted in Fig. 5. We also computed the PDFs ex¬ 
pected for pure shape noise data using 1000 random re¬ 
alizations of “noise K. maps”, and the mean and RMS of 
1000 PDFs are also plotted in the figure. The upper panel 
of Fig. 5 shows the peak PDF measured from /C maps, 
where clear deviations from pure shape noise PDF are 
seen. The deviations originate mainly from three sources: 
( 1 ) massive clusters of galaxies that account for very high 
peaks (say S/N > 4.5), (Miyazaki et al. 2002a; Hamana 
et al. 2004) which we shall go into details later, (2) cos¬ 
mic large-scale structures which contribute to peaks at 
5'/iV^ 4.5 (Maturi et al. 2010; Marian et al. 2010), and (3) 
residual systematics. In order to estimate the impact of 
the last one, we use the B-mode 1C map (/C_b), which is the 
same smoothed convergence map as /C but obtained from 
45-degree rotated shear. In the standard gravity theory, 
only a small amount of B-mode lensing signal is produced 
(Schneider et al. 2002; Hilbert et al. 2009) which should be 
negligible for the current analysis. Since in the absence of 
lensing-induced B-mode signal, the /Cs-peak PDF should 
follow the pure shape noise PDF, it can be used for a diag¬ 
nosis of residual systematics. It is found from the bottom 
panel of Fig. 5 that the /C^-peak PDF deviates from the 
pure shape noise PDF, indicating existence of some resid¬ 
ual systematics. Although a cause of the excess /Cs-peaks 
is unknown, it should be safer to consider that a simi¬ 
lar amount of non-lensing-originated peaks may exist in 
the /C maps. It can be also said from the ICb peak PDF 
that the probability of having a non-lensing-induced peak 
with S/N > 5 is very rare. In addition to this fact, with a 
limited number of peaks with S/N > 5, it is possible to in¬ 
spect all the individual high peaks. Therefore we decided 
to use only the peaks with S/N > 5 for our cosmological 
study in the next section. Understanding the origin of the 
excess ICb peak PDF must be a very important subject 
for future studies of weak lensing cluster counts using data 
from large surveys, however it is beyond the scope of this 
paper and we leave it for a future work. 

3.2. High S/N peaks 

/C-peaks with S/N > 4.5 are marked with open cir¬ 
cles in Fig. 4, and are summarized in Table 2. We 
searched known cluster database taken from a compila¬ 
tion by NASA/IPAC Extragalactic Database (NED)^ for 
possible luminous baryonic counterparts within 2 arcmin 
from peak positions, and matched clusters of galaxies are 
also presented in Table 2. Most of the matched known 
clusters were identified by galaxy concentration in opti¬ 
cal/infrared multi-band data and photometric redshifts of 
the clusters were evaluated. For some cases the spectro¬ 
scopic redshifts of the brightest cluster galaxy (BCG) were 
measured. It turns out that all but XMM-LSS I, which 
has the peak height of S/N = 4.7 and is located in “edge- 
region” , have matched known clusters. Although for some 

^ http://ned.ipac.caltech.edu/ 



^shape 


Fig. 6. The probability distribution of the shape noise (Tshape 

evaluated evaluated on each grid point using 1000 realiza¬ 
tions of “noise /C maps” as described in §3.1. The mean and 
RMS of this distribution are found to be (cTghape) = 0.017 and 
RMS(o-gi,ape) = 0.0016, respectively. 

peaks, matched clusters are located at slightly discrepant 
positions {9sep > 1 arcmin), all the peaks with S/N > 5 lo¬ 
cated in the “data-region” have a relatively well matched 
cluster of galaxies (0sep < 1 arcmin). Photometric redshifts 
of those closely matched counterparts are within the high 
lensing efficiency range (0.1 < 2 < 0.5), supporting an ac¬ 
tual correspondence between galaxy over-density and high 
weak lensing signal. We thus conclude that the sample of 
S/N > 5 peaks located in the “data-region” are not con¬ 
taminated by a false signal, and we use the number counts 
of Npeak{S/N > 5) = 6 in an effective area (“data-region” 
only) of 9.0 degree^ for the cosmological study in the next 
section. 

4. Constraints on the cosmological parameters 
and M — c relation 

Here we place constraints on the cosmological parame¬ 
ters (erg and Clm) and the mass-concentration {M — c) re¬ 
lation using the measured number counts of weak lensing 
clusters. To this end, we compare the measured counts 
with the theoretical model in the standard method. 
Below, we first describe the theoretical model of weak 
lensing cluster counts, then next we evaluate the error in 
the measured number counts which consist of the Poisson 
error and the cosmic variance, and finally we present re¬ 
sults. 

4 . 1 . Theoretical model of weak lensing peak counts 

We adopt the theoretical model of the weak lensing peak 
counts developed in Hamana et al. (2012), which is based 
on halo models. We refer the reader to the above reference 
for details of model and its comparison with mock numer¬ 
ical simulations. In short, the model is constructed under 
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Table 2. Summary of high peaks with S/N > 4.5. 


Field - ID 

R.A. 

[deg] 

Dec. 

[deg] 

Peak 

S/N 

flag® 

Q h 

f'sep 

Matched known clusters (0sep < 2 aremin) 

Zs,bcg'^ Name and reference 

XMM-LSS 1 

36.6087 

-5.0671 

4.7 

edge 





GOSMOS 1 

150.3017 

1.5838 

4.8 

- 

1 .2' 

0.36 


#59 of (1) 






1.5' 

0.242 


NSG J100113-f013335(2) 






1.5' 

0.37 


LSS 27 of (3) 






1.9' 

0.3610 

0.3638 

GOSMOS GL J100112.4-H013401W 

GOSMOS 2 

149.4210 

2.5585 

5.1 

edge 

0.5' 

0.3787 


WHL J095742.9-f023332(5) 






1 .0' 

0.17 


LSS 36 of (3) 






1.4' 

0.362 

0.3733 

GMBGG J149.40425-f02.5738l(®) 






1.4' 

0.39 


LSS 26 of (3) 

Lockman-hole 1 

161.6291 

59.6266 

7.6 

edge 

0.7' 

0.2317 

0.2282 

WHL J104625.5-f593736(5) 






1 .0' 

0.5672 


SWIRE GL J104635.5-f593732(4) 

ELAIS-Nl 1 

244.2695 

53.6936 

4.6 

- 

1 .8' 

0.9061 


SWIRE GL J161710.9-f534217('‘) 

ELAIS-Nl 2 

242.8739 

54.2766 

5.2 

- 

0.5' 

0.3794 


NSG 4161132-^541649(2) 






0 .8' 



SWIRE GL 4161130.1-^541704(4) 






0.9' 

0.3307 


WHL 4161135.9-^541634(5) 

ELAIS-Nl 3 

242.8169 

55.5667 

6.7 

- 

0 .8' 

0.2246 


SWIRE GL 4161118.8-f55343l(4) 






1 .1' 

0.2156 


WHL 4161121.1-^553308(5) 






1 .8' 



ZwGl 1610.0-f5543(^) 

ELAIS-Nl 4 

242.6503 

54.0947 

5.1 

- 

1 .0' 

0.3786 

0.3370 

SWIRE GL 4161045.2-^540635(4) 






1 .2' 

0.3618 


2XMM 4161040.5-f540638(®) 

ELAIS-Nl 5 

241.6525 

53.8947 

5.3 

- 

0 .2' 

0.3432 

0.3505 

SWIRE GL 4161045.2-^540635(4) 






0.7' 

0.348 


GMBGG 4241.67338-^53.89618(5) 

ELAIS-Nl 6 

241.2632 

55.6340 

5.3 

- 

0 .1' 

0.25385 


MaxBGG 4241.26104-^55.63254(9) 






1 .0' 

0.275 


WHL 4160503.4-^553703(5) 






1 .2' 

0.2299 


NSG 4160459-f553700(2) 

ELAIS-Nl 7 

241.1661 

54.5305 

8.2 

- 

0 .2' 

0.2483 


WHL 4160439.4-^543136(5) 


“ Flag for edge regions (see §3.1 for the definition of the edge region). 

^ Separation between the K. peak and the cluster position given in the references. 

Photometric redshift of clusters given in the references. 

^ Spectroscopic redshift of BCG cluster given in the references. 

Reference list: ^^^^Finoguenov et al. (2007), ^^^Gal et al. (2003), ^^^Scoville et al. (2007), ^^^Wen & Han (2011), ^^^Wen 
et al. (2009), ^^^Hao et al. (2010), ^^^Zwicky et al. (1961), ^^^Takey et al. (2011), ^^^Koester et al. (2007). 


the assumption that high peaks are dominated by lensing 
signals from single massive haloes, for which it is assumed 
that the density profile is described by the truncated NFW 
model (Navarro et al. 1997; Baltz et al. 2009) and mass 
function is given by the model of Warren et al. (2006). 
Effects of intrinsic galaxy shape, the diversity of dark mat¬ 
ter distributions within haloes, and large-scale structures 
are included in approximate approaches which were tested 
against the mock numerical simulations. Model parame¬ 
ters that should be specified for making the theoretical 
prediction of the weak lensing cluster counts are the cos¬ 
mological parameters, M — c relation of the dark matter 
haloes, the shape noise CTshape, and the redshift distribu¬ 
tion of source galaxies. Below we will describe them in 
turn. 

We treat and erg as free parameters, and fix the 
other cosmological parameters as follows: Hubble param¬ 
eter h = 0.7, the cosmological constant Ha = 1 — the 
baryon density Hf, = 0.051, and the spectral index n = 0.96, 
which are all consistent with the recent GMB experiments 


(Hinshaw et al. 2013; Planck Gollaboration et al. 2014). 
For the M — c relation, we assume the following form, 

/ Jut . \ -0 075 

= (1 + .)-°^ (3) 

which is motivated by the hndings from dark matter N- 
body simulations by Duffy et al. (2008); Maccio et al. 
(2008); Klypin et al. (2011); Prada et al. (2012), but we 
treat the normalization as a free parameter. Klypin et al. 
(2011) found Co = 9.6 for the cosmological model with 
Hm = 0.27 and erg = 0.82. However, it should be noted 
that the mass-concentration relationship is known to be 
dependent on the cosmological model (Maccio et al. 2008; 
Zhao et al. 2009; Kwan et al. 2013; Diemer & Kravtsov 
2015). 

The peak fC value from a halo is computed by assuming 
the tangential shear profile of a truncated NFW halo con¬ 
volved with the smoothing function, eq. (2). Then it is 
converted to S/N by specifying the shape noise, Ushape- 
Since Ushape depends on the number density of source 
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Fig. 7. Redshift distribution of galaxies used for the weak 
lensing mass reconstruction derived adopting the COSMOS 
photometric redshift catalog (Ilbert et al. 2009). The solid 
curve shows the fitting model by equation (8) with the pa¬ 
rameters denoted in the panel. 

galaxies, it fluctuates depending on position. As we have 
described in §3.1, we evaluated ctshape on each grid of weak 
lensing mass maps. We computed the probability distribu¬ 
tion of the shape noise, which is presented in Fig. 6. The 
distribution is not very widely spread: its mean and RMS 
are found to be (ashape) = 0.017 and RMS(CTshape) = 0.0016, 
respectively. We take this effect into account in the fol¬ 
lowing way, 

peak (*5*/ JV') = J' d(7shapeP{^shape')^peak{S /fV 5 (Tshape) ; (4) 

where iVpeak(<S'/A^; o-ghape) is the theoretical prediction of 
the weak lensing cluster number counts for a fixed fTghape 
and P(crshape) is the probability distribution function. 
Note that crghape measured from the data is affected by 
the multiplicative bias in shear measurement (§2.2). In 
order to take it into account, we calibrate the shape noise 
RMS by CTshape = CTshape/ (1 - {m)), and P(CTshape) with the 
measured distribution function shown in Fig. 6 is used to 
evaluate eq. 4. 

In order to infer the redshift distribution of source 
galaxies, we utilized the COSMOS held data. We merged 
the source galaxy catalog used for the weak lensing mass 
reconstruction with the public COSMOS photometric red¬ 
shift catalog (Ilbert et al. 2009), and computed the red¬ 
shift distribution by adopting the photometric redshift. 
The derived redshift distribution is presented in Fig. 7. 
We ht the redshift distribution with the following func¬ 
tion (Fu et al. 2008): 

-a _|_ 2'^b 

= ^ + 

where the normalization, was determined by imposing 


Fig. 8. The probability distribution function of the number 
counts of peaks with S/N > 5 found in 9.0 degree^ fields. 
The distribution is derived from 10,000 realizations randomly 
taken from the full-sky mock weak lensing data. The mean 
and RMS of the distribution is 7.03 and 3.26, respectively. 

/dz ns{z) = 1 within the integration range 0 < 2 < 6. We 
found that a parameter set (a,6,c) = (0.4,5.5,1.5) gives a 
reasonably good fit to the data, as shown in Fig. 7. The 
mean redshift of the model is 1.10 which is in very good 
agreement with the value derived from the photometric 
redshift data (1.09). In order to examine the variation in 
the redshift distribution due to the difference in the expo¬ 
sure time, we divide the COSMOS field into two subfields: 
one is the central region (^0.5 degree^) where the expo¬ 
sure time is longer than 3,600 sec (with a median of 4,920 
sec), and the other is the outer region where the exposure 
time is shorter than 3,600 sec (but mostly 2,400 sec) (see 
Taniguchi et al. 2007, for variation of the limiting mag¬ 
nitude). We computed the redshift distributions of those 
two samples separately, and found that the distributions 
are nearly indistinguishable, with mean redshifts of 1.093 
and 1.090 for the inner-deep and outer-shallow samples, 
respectively. This may be partly due to our selection crite¬ 
ria that i' < 25 [AB mag], detection S/N > 8 and FWHM 
larger than the stellar FWHM, which may result in rejec¬ 
tion of small and faint distant galaxies. According to this 
finding and the fact that the redshift information for the 
other three fields is not available, we do not consider a 
possible field-to-field variation in the redshift distribution 
but use eq. (5) for computation of the theoretical predic¬ 
tion. 

4 . 2 . Uncertainty in the number counts 

In order to evaluate the uncertainty in the weak lensing 
cluster counts, we utilize mock weak lensing data gener¬ 
ated from full-sky gravitational lensing ray-tracing simu¬ 
lations. Details of the simulation are given in Shirasaki, 
Hamana & Yoshida (2015 in prep.), here we describe 
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only aspects that are directly relevant for this work. 
We adopted the cosmological model with ftm = 0.278, 
fiA = 0.722, (Tg = 0.82, h = 0.70 and n = 0.97. These 
parameters are consistent with the WMAP nine-year re¬ 
sults (Hinshaw et al. 2013). Gravitational leasing ray¬ 
tracing was performed through the dark matter distri¬ 
bution generated by the Wbody simulations. In order 
to cover full-sky up to 2 > 1, we combined, in a nested 
structure, 8 simulation boxes of different sizes with box 
side length from 450h“^Mpc to 3,600h“^Mpc. The N- 
body simulations were run with 2048^ dark matter par¬ 
ticles using the parallel Tree-Particle Mesh code Gadget2 
(Springel 2005). Light ray paths and the magnification 
matrix were evaluated by the multiple-plane ray-tracing 
algorithm of Teyssier et al. (2009) and Becker (2013). 
We used HEALPix pixelization of a sphere (Gorski et al. 
2005). The angular resolution parameter nside was set 
to 8192. The corresponding angular resolution is 0.43 ar- 
cmin. We adopted a fixed source redshift of Zs = 1.09. 
On each pixel center, we computed the lensing shear. 
We added a random galaxy shape noise drawn from the 
truncated two-dimensional Gaussian to the lensing shear 
data (Hamana et al. 2012). The a of the Gaussian func¬ 
tion was determined so that the shape noise in the resul¬ 
tant /C map became equal to the mean observed value 
(after calibration with the multiplicative bias), that is 
dshape = Cshape/(1 — (w)) = 0.018. We adopted the same 
smoothing function as the actual measurement, eq. (2), 
and computed K, from the noise-added shear data. Since 
we used the dark matter distribution from cosmological 
A^-body simulations, the influence of large-scale struc¬ 
tures on weak lensing cluster counts should be properly 
included. 

Having generated the mock fC data, we count the num¬ 
ber of peaks with the threshold of S/N > 5 in a contigu¬ 
ous 9.0 degree^ region, which is randomly taken from the 
full-sky simulation data. We repeat this counting 10,000 
times, allowing overlap in regions between different real¬ 
izations. The probability distribution function of the re¬ 
sulting number counts is shown in Fig. 8. The mean and 
RMS of the distribution was found to be 7.03 and 3.26, 
respectively. Note that the solid curve in Fig. 8 represents 
the Gaussian distribution with the same mean and RMS, 
demonstrating that the measured distribution can be rea¬ 
sonably approximated by the Gaussian distribution. We 
may safely assume that the RMS consists of Poisson error 
and cosmic variance. Taking into account the difference 
between the observed number counts of 6 and the mean 
number counts of 7.03, we adopt ctat = 3.1 for the error 
in the observed weak lensing cluster counts, and thus we 
finally have A^peak = 6 ± 3.1. 

4-3. Results 

In order to place constraints on cosmological parameters 
and the M — c relation using weak lensing cluster counts, 
we compare the measured counts with the theoretical pre¬ 
diction utilizing the standard method. Note that the 
resulting x^ value does not represent the goodness-of-fit, 
as we have only one observed value, Npeak = 6±3.1. Thus 


Fig. 9. Constraints on ^rn erg derived from the weak 
lensing cluster counts. Here the normalization of the M — c 
relation is assumed to be cq = 9.6 (Klypin et al. 2011). Gray 
scale shows the with the inner green and outer red (up¬ 
per limit only) contours indicating 68.3% (x^ =2.30), 95.4% 
(^2 _ g confidence levels, respectively. 



0.6 0.8 1 1.2 


Fig. 10. Same as Fig. 9 but on the erg and cq. Here the 

density parameter was fixed with flrn = 0.278. 

we use the x^ value only for placing a confidence region 
on parameter spaces, adopting confidence levels assuming 
a normal distribution. We treat ug and cq as free 
parameters. Constraints obtained from the weak lensing 
cluster counts have degeneracy between those parameters 
(Mainini & Romano 2014; Cardone et al. 2014). However, 
since we have only one observed value with poor statis¬ 
tics, we do not explore the three-parameter space, but we 
examine two-parameter space, fixing one remaining pa¬ 
rameter. 

Figure 9 shows the derived constraint on — ug plane. 
Here we did not treat co as a fitting parameter but adopted 
the value (cq = 9.6) found from dark matter V-body simu¬ 
lations by Klypin et al. (2011) in which the ACDM model 
with ilm = 0.27 and ug = 0.82 was employed. Therefore 
the constraint mainly comes from the cosmology depen¬ 
dence of the dark matter halo mass function with a minor 
contribution from the cosmology dependent lensing effi¬ 
ciency. Accordingly, the constraint exhibits the ilm — <7g 
degeneracy in a similar fashion to X-ray cluster counts 
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Fig. 11. Constraint on a single parameter cq derived from 
the weak lensing cluster counts. Here Qrn and as are fixed 
to be the best fitting model of either the Planck (solid curve 
Planck Collaboration et al. 2014) or WMAP 9-year (dashed 
curve; Hinshaw et ah 2013). 

(Borgani et al. 2001; Henry 2004; Vikhlinin et al. 2009; 
Bohringer et al. 2014). 

Figure 10 shows the derived constraint on the as — cq 
plane, where the density parameter was fixed with flm = 
0.278. The degeneracy in the as — cq space results from 
two relationships: one is increasing dark matter halo num¬ 
ber density for larger erg, and the other is higher 5'/7V-peak 
for larger cq. Both of them result in larger weak lensing 
cluster counts. The upper limit on cq is hardly placed, 
because of the degeneracy in erg — Cq and, for a given 
smoothing kernel, an excessively concentrated halo com¬ 
pared with the smoothing scale causes little variation in 
the peak height (Hamana et al. 2004). Although the cur¬ 
rent constraint does not even place a useful upper/lower 
limit on those parameters due to the poor statistics, this 
demonstrates that weak lensing cluster counts can be a 
useful probe of the cosmic structure formation scenario. 
The same argument, but based on a slightly different 
methodology, was already made by Liu et al. (2014b). 

Finally, for a demonstrative purpose, we present in 
Figure 11, the constraint on a single parameter cq for two 
representative cosmological models: the Planck = 
0.315, as = 0.83) (Planck Collaboration et al. 2014) and 
WMAP 9-year {flm = 0.278, erg = 0.82) (Hinshaw et al. 
2013). It is obvious that the measured weak lensing cluster 
counts are consistent with the currently popular ACDM 
model. Specifically, the result from V-body simulations 
by Klypin et al. (2011) (cq = 9.6) is well within the 68.3% 
confidence region of both the Planck and WMAP9 cos¬ 
mologies. 


5. Summary and discussions 

We presented the results of weak lensing cluster counts 
obtained from 11 degree^ Subaru/SuprimeCam data. 
For the first time, we used the weak lensing cluster 
counts to place cosmological parameters and the mass- 
concentration relationship. In doing this, we explored, in 
an experimental manner, the capability and current issues 
of weak lensing cluster counts for cosmological studies. 

The data were selected under the conditions that the 
exposure time longer than 1,800 sec along with the seeing 
condition of FWHM < 0.7 arcsec. These selection criteria 
result in the effective number density of galaxies usable 
for weak lensing shape measurement of Ug ^ 20 arernin”^, 
which is about twice as large as those of previous works 
dealing with weak lensing peak statistics, namely, the 
CFHTLenS by Liu et al. (2014a) and CFHT Stripe-82 
Survey by Liu et al. (2014b). Since the shape noise in the 
1C map scales as ae oc Ug higher galaxy density lifts 
up the peak height of a lensing signal of massive clusters. 
This allows us to adopt a high threshold S/N, in our case 
S/N > 5. With high S/N, the contamination rate due to 
false signals is small, which enables us to construct a weak 
lensing-selected cluster sample. 

We used lens fit for weak lensing shear measurement. 
We examined the additive bias in the measured galaxy 
ellipticities by evaluating weighted mean ellipticities as a 
function of galaxy size and magnitude. We found non¬ 
zero additive bias for both ei and 62 , which were found 
to be mainly caused by faint and small galaxies, as shown 
in Fig. 3. This is in marked contrast to the result of 
CFHTLenS in which small and bright galaxies were found 
to be a cause of additive bias in 62 (Heymans et al. 2012). 
The cause of these non-zero biases is unclear, but those 
different trends for different data sets may suggest that it 
may be caused by some instrumental or data-dependent 
effect. We evaluated the amplitude of the multiplicative 
bias, adopting the empirical fitting function derived from 
CFHTLenS image simulations (Miller et al. 2013), and 
found the weighted average value of (m) = —0.059. This 
does not have influence on the measured peak-S'/V, as 
both the signal and noise RMS are affected in the same 
way, and thus those are canceled out. However, this must 
be properly taken into account in making the theoretical 
prediction. Since in the current study we considered the 
statistical quantity (peak counts), we adopted a statistical 
correction scheme that calibrating the shape noise RMS 
with the average bias factor (1 — (m)). 

We found 6 peaks with S/N > 5. For all the peaks, pre¬ 
viously identified clusters of galaxies are matched within 
a separation of 6 sep < I aremin. Most of the matched 
clusters were identified as a galaxy concentration in the 
sky and in multi-band color space of optical/infrared data. 
This provides us with a firm inspection of the purity of the 
weak lensing cluster sample. In fact, optical multi-band 
images are well suited to cross-identification, as the sensi¬ 
tivity of weak lensing mass map to cluster finding is most 
effective for clusters within the redshift range 0.1 ^ z ^ 0 . 6 , 
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where the red-sequence of cluster galaxies is located in op¬ 
tical color space. Thus, ongoing optical multi-band galaxy 
surveys with a single band observation being optimized 
for weak leasing shape measurement, namely the Kilo- 
Degree Survey (KiDS)^, Dark Energy Survey (DES)"^, and 
Hyper-SuprimeCam (HSC) Survey (Miyazaki et al. 2012), 
may provide self-contained datasets for weak leasing clus¬ 
ter finding. Also, photometric redshifts will be estimated 
from multi-band data for each galaxy, which enable the 
source galaxy redshift distributions to be accurately de¬ 
termined. 

We evaluated the statistical error in the weak leas¬ 
ing cluster counts using mock weak leasing data gen¬ 
erated from full-sky ray-tracing simulations, and found 
A^peak = 6 ± 3.1. This error consists of the Poisson noise 
and the cosmic variance, with a larger contribution from 
the former. 

We compared the measured weak leasing cluster counts 
with the theoretical model developed by Hamana et al. 
(2004, 2012), in which the effects of intrinsic galaxy 
shape, the diversity of dark matter distributions within 
haloes, and large-scale structures are included. Utilizing 
the standard method, we placed constraints on the 
~ erg plane, which was found to be consistent with 
currently standard ACDM models such as Planck model 
(Planck Collaboration et al. 2014) and WMAP 9-year 
model (Hinshaw et al. 2013), though the constraints are 
much broader than those CMB experiments due to the 
poor statistics (A^peak = 6 ±3.1). It was demonstrated 
that the weak leasing cluster counts may place a unique 
constraint on ug — cq plane. 

Einally we discuss prospects for ongoing/future surveys, 
taking the DES and HSC surveys as examples. Even 
taking into account a lower usable galaxy number den¬ 
sity for those surveys, due to the shallower survey depth 
than this study, it may be feasible to obtain Apeak ^ 400 
clusters with S/N ^ 5 (which can be achieved with ~ 0.3 
peak/degree^ for HSC survey of 1,400 degree^, and ^ 0.08 
peak/degree^ for DES of 5,000 degree^). In that case, 
the error is dominated by the Poisson error and the ex- 
pected fractional error will be ^ 0.05. This is 10 

times better than the current study, 3.1/6 ^ 0.5. Therefore 
confidence regions become smaller and thus tighter con¬ 
straints will be obtained; taking Eigure 11 as an illus¬ 
trative example, the corresponding 68.3% confidence level 
becomes = 0.1. This demonstrates that this method¬ 
ology will be a useful probe of the cosmology and M — c 
relation, though, in reality, marginalization over other pa¬ 
rameters and degeneracy among parameters weaken con¬ 
straints. More quantitative investigations of the capability 
of this method based on a Eisher analysis were made by 
Cardone et al. (2014); Mainini & Romano (2014). Having 
found that tight cosmological constraints can be achieved 
using the significantly improved statistics expected for on¬ 
going/future surveys, it is necessary to check whether the 
theoretical model is correspondingly accurate, which is be- 

® http://kids.strw.leidenuniv.nl 
http: //www. darkenergysurvey. org 
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yond the scope of this paper and is left for a future work. 
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Appendix 1. Method of conversion from SIP to 
TPV 

As described in § 2 , the instrumental distortion of im¬ 
ages was corrected in the process of mosaic stacking along 
with astrometric calibration of CCD data using hscPipe 
in which the relationship between pixel coordinates and 
world coordinates systems (WCSs) is represented by SIP 
convention (Shupe et al. 2005). In lensfit, on the other 
hand, TPV convention® is implemented. Here we describe 
the method for conversion of WCSs from SIP to TPV. 

A. 1.1. SIP convention 

First, we summarize the SIP convention following Shupe 
et al. (2005). Let u^v be relative pixel coordinates on a 
detector device with origin at CRPIXl, CRPIX2, and x',y' 
be distortion corrected “intermediate world coordinates” 
in degree with origin at CRVALl, CRVAL2. Then 

/ CD1_1 cmj2\( u + f{u,v)\ , , 

CD2_1 CD2.2 y v + g{u,v) 

where CDz_j matrix encodes scaling, rotation and skew, 
and f{u,v) and g{u,v) are polynomial functions that rep¬ 
resent distortion, given by 


f{u,v) = ^A.Aj u^v\ 

i^j 

* + J < A.ORDER, 

(A2) 

g{u,v) = ^^B.i.j u^v^, 

* + J < B.ORDER. 

(A3) 


ij 


In the standard SIP convention, A_0RDER and B_0RDER are 
allowed to range from 2 to 9. 

A. 1.2. TPV convention 



A.1.3. Conversion from SIP to TPV 


We derive conversion equations from SIP coefficients to 
TPV coefficients. To do so, we limit the polynomial order 
of SIP (A_0RDER and B_0RDER) being equal to or less than 
7th, because that of TPV is defined up to 7th. 

We can set PVCj coefficients for r"-terms® zero, be¬ 
cause there are no corresponding terms in SIP polynomi¬ 
als. Also, we can safely use the same origins of {u,v) and 
(x',t/') coordinates for TPV as those for SIP, because we 
derive the same distortion polynomial models as those of 
SIP but in the different (TPV) convention (thus the same 
CRPIXl, CRPIX2, CRVALl and CRVAL2 can be used). Then, 
we can set PV1_0 = PV2_0 = 0 as there are no corre¬ 
sponding terms in SIP. In addition, utilizing redundant 
degree of freedom in TPV convention, we can safely take 
the same CDi_j matrix as one of SIP. 

Equating the SIP and TPV coefficients of each 
terms by order-by-order, we have the following simultane¬ 
ous equations: For the 1st order, 


CD1_1 

CD1_2 

CD2_1 

CD2_2 


CD 1.1 
CD 1.2 

CD2.1 

CD2.2 


CD2.1 

CD2.2 

CDl.l 

CD1.2 


PVl.l 

PV1.2 

PV2.1 

PV2.2 


,(A7) 

,(A 8 ) 


which can be inverted and can be solved as, 


f PVl.l \ _ f PV2.1 \ _ f 1 \ 
PV1.2 ) ~ \ PV2.2 ) ~ \ 0 )' 


(A9) 


These are a natural consequence of distortion polynomial 
functions of SIP convention being defined only for 2nd and 
higher order terms (see eqs. (Al) and (A2)). For the 2nd 
order, we have 


Next, we summarize the TPV convention, in which rel¬ 
ative pixel coordinates (u, v) are first converted to “dis¬ 
torted world coordinates” {x,y) by 

fx\ /CDl.l CD1.2 \ / u \ 

( 9 j = ( CD2.1 CD2J j ( » j ■ 

Then “intermediate world coordinates” {x',y') can be 
given by polynomial functions of x and y as, 

x' = PVl.O 4- PVl.l X + PV1.2 y + PV1.3 r 
+PV1.4 x^ + PV1.5 xy + PV1.6 
-PPVl.7 a;® -P PV1.8 x'^y + PV1.9 xy"^ 

-PPVl.lO 2/® + PVl.ll (A5) 

y' = PV2.0 -P PV2.1 y -P PV2.2 x + PV2.3 r 
+PV2.4 y^ + PV2.5 xy + PV2.6 x^ 

-PPV2.7 -P PV2.8 xy^ + PV2.9 x^y 

+PV2.10 X® + PV2.11 r® -P • • •, (A6) 

where r = y/z?-P?/^, and the TPV polynomial functions 
have odd-power terms of r. The polynomial order of TPV 
is defined up to 7th, thus there are 2 x 40 PVi.j coefficients 
including 2x4 odd r^-terms. 


/ A.2.0 CDl.l-PB.2.0 CD1.2 \ 

A.1.1 CDl.l-PB.1.1 CD1.2 = 

\ A.0.2 CDl.l-PB.0.2 CD1.2 j 

/ CD1.i2 CD1.1CD2.1 

2CD1.1CD1.2 CD1.1CD2.2-I-CD1.2CD2.1 

\ CDI.22 CD1.2CD2.2 

/ PV1.4 \ 

X PV1.5 

\ PV1.6 j 

and 

/ A.2.0 CD2.1 -P B.2.0 CD2.2 \ 

A.1.1 CD2.1 -P B.1.1 CD2.2 = 

\ A.0.2 CD2.1 -P B.0.2 CD2.2 ) 

/ CD2.12 CD2.1CD1.1 

2CD2.1CD2.2 CD2.1CD1.2-I-CD2.2CD1.1 

y CD2.22 CD2.2CD1.2 

/ PV2.4 \ 

X PV2.5 

\ PV2.6 j 

In the same manner, one can derive simultaneous equa¬ 
tions in matrix form for higher order terms. Therefore, 


CD2.D 

2CD2.1CD2.2 

002 . 2 ^ 


(AlO) 


CDl.D 

2CD1.1CD1.2 

001 . 2 = 


(All) 


http: //fits.gsfc.nasa.gov/registry/tpvwcs/tpv.html 


® Those are PV1_3 , PV1_8 , PV1_23 , PV1_39 , PV2_3 , PV2_8 , 
PV2_23 and PV2_39 . 
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TPV coefficients, PYi-j , can be obtained by solving those 
simultaneous equations using the standard matrix inver¬ 
sion method. In an actual computation, we performed 
matrix inversion and multiplication numerically. 


